Geometrically Frustrated Crystals: Elastic Theory and Dislocations 
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Elastic theory of ring-(or cylinder-) shaped crystals is constructed and the generation of edge 
dislocations due to geometrical frustration caused by the bending is studied. The analogy to su- 
perconducting (or superfluid) vortex state is pointed out and the phase diagram of the ring-crystal, 
which depends on radius and thickness, is discussed. 
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What is the lowest energy state of matter? The an- 
swer to this fundamental question has been believed to 
be the "crystalline state" as far as the system has an 
infinite extension and quantum fluctuation is negligible. 
Most elements or compounds freeze into crystals when 
cooled down slowly. Even the electrons form the so- 
called Wigner crystal if the density is low enough. The 
basic character of the crystalline state is the breaking of 
translational symmetry; the system is divided into small 
fundamental cells. Although some exceptions, such as 
Penrose lattice or quasicrystals, have been known, this 
profound form of existence occupies the most intriguing 
expositions in the "mineralogical department of natural 
history museum" . 

In recent decades, with the help of developed fine pro- 
cessing technology, especially symbolized by the term 
"nanotechnology" , a new field seems to be added to 
the studies of crystallography. The topic of this field 
is the crystallization in nano-scales. A group of mate- 
rials which is most attracting attention from both basic 
and applied scientists is evidently Fullerenes and nan- 
otubes 0, 0, S IE IE 0) where spherical or cylindrical 
form of crystals and their subspecies are found to be 
stable structures of carbon in nano-scales. Tanda et al. 
have shown that quasi-one-dimensional crystals, such as 
TaS3 or NbSe3 , which are usually grown as needle-shaped 
crystals, can be made into rings or cylinders by control- 
ling the growing conditions 

[IE EE III 03 They have 
even succeeded in producing single-crystal Mobius strips. 
These crystals with topologically nontrivial shapes are 
named "topological matters". Some other cylindrical 
crystals were also found recently, such as M0S2 ,13] . The 
expanding variety of the topological matters seems to be 
invoking questions about the nature of crystallization in 
nano-scales. 

The crystal structure largely affects the physical prop- 
erties of the system and, therefore, these crystals pro- 
vide excellent experimental environments in investigating 
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topological physics. For example, TaS3 and NbSe3 show 
superconductivity and charge density wave order under 
appropriate conditions [l4|, arK1 , then, topology-induced 
novel states are expected. Several theoretical works have 
been elaborated in this direction [TEl Ibrj ITtI Il8| . and 
experimental investigations are under way. 
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FIG. 1: The structure of the ring-(or cylinder-)shaped crystal 
made of a orthorhombic lattice. Because of the expansion and 
contraction caused by the bending, edge dislocations are likely 
to be generated. Here W is the thickness of the ring and R is 
the radius. R is defined at the center of the thickness. 

In order to understand the physical properties of the 
"topological matters" , the knowledge about crystal struc- 
ture is indispensable. In this paper we study the struc- 
tures of ring-(or cylinder-) shaped crystals, which we call 
simply the "ring crystals" hereafter. Special attention 
is paid to the formation of crystal defects due to the 
elastic stress which topological matters undergo because 
of their nontrivial forms. In case of the ring crystals the 
edge dislocations are likely to be generated because of the 
bending, as one can see from Fig. In this paper we 
derive, within the linear elastic theory, the condition for 
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the formation of edge dislocations in a ring crystal whose 
radius and thickness are R and W, respectively. Analogy 
between the superconducting (or superfluid) vortices and 
the dislocations in a ring crystal is invoked. Our results 
may be useful in understanding physical phenomena in 
ring crystals. 

Elastic theory of a ring crystal 

To construct the elastic theory of a ring crystal we 
start by making a rectangular strip of a crystal into a 
ring as depicted in Fig. [3 which introduces elastic de- 
formation. We assume that the system is uniform in 
z-direction (along the axis of the cylinder) and suppress 
the z-coordinatc hereafter. The free energy is considered 
to be that for unit length in z-direction. The crystal is 
anisotropic and it is assumed that the elastic coupling in 
the direction along A-C and B-D is much stronger than 
that along A-B and C-D (see Fig. EJ. This situation is 
realized in most ring crystals grown until now. 



-1/2 







„ - 

/ 





L/2 
— *► 



X 




FIG. 2: Making a rectangular strip of a crystal into a ring. 
The strip (a) is made into a ring (b) by introducing the elastic 
deformation. The system is assumed to be uniform in z-axis, 
which is perpendicular to the page. 

The elastic free energy of the ring crystal is constructed 
by examining how the linear element of the original crys- 
tal df = [dx, dy) is modified by the bending in the ring 
crystal. Let us suppose that the position in the original 
crystal r = (x, y) is moved to R = (X, Y) in the ring 
crystal. Here orthogonal coordinate system is used. We 
obtain the following relation, 

,„ Trt f, . x + u , . x + u\ 

(X,Y) = I (r + v)cos— —,(r + v)sm— — J , (1) 

where r = r(y) denotes the radius of the lattice plane 
whose y-coordinate is y in the original system, and 
u = u(x,y) and v = v(x,y) are small deviations around 
the complete ring shape (i.e., dislocation-free ring). We 



define the metric tensor by (dR) 2 = J^i j gijdxidxj, 

where g tj = §£-■§£- = {1, 2}, and xi = x,x 2 = y). 
The elements of the metric tensor are given as 
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912 = 9x V = -jp(r + v) 2 {l + u x )u y + v x (r' + v y ), (4) 

where r' — dr/dy, u x = du/dx ■ ■ ■ etc. In this paper, we 
assume r(y) = R + y (-W/2 < y < W/2) for simplicity. 
Strain tensor Uij, which describes the elastic deforma- 
tion of the crystal, is given by = \{gij — Sij), where 
Sij is the Kronecker's delta. Assuming an orthorhombic 
system for simplicity, the elastic free energy is given by 
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where L = 2nR is the circumference of the ring [lflj . 
We also assume that the principal axes of the crystal 
are along x- and j/-axis, which yields \ xxyy — —X xyxy /2. 



From now on we denote Ai = X X xxx : A2 
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A3 — X x y X y X xx yy /2. 

Next we determine the free energy of the ring crystal 
within the linear elasticity. To do this we expand the free 
energy with respect to small parameters, u x , u y , v x , v y , 
y/R and v/R to the second order, which yields 
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Here we limit ourselves to the crystals with W <C R- We 
further neglect w-field here, since A2 3> A3 in most ring 
crystals, and dominant elastic fluctuation is expected to 
arise not from v but from u. Then we obtain the following 
free energy, 
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where d x and cf>Q are the lattice constant in ^-direction 
and the superconducting magnetic flux quantum, respec- 
tively. This is nothing but the free energy of an isotropic 
superconductor poj . The superconducting "phase" 9 un- 
dergoes the quantization condition J c dl ■ V# = 2-jm v , 
where C is a simple loop and n v is the total number of 
vortices enclosed in C . The vortices correspond to the 
edge dislocations in the crystal system. Thus the quan- 
tization of vorticity corresponds to the quantization of 
the Burgers vector. Here the vector potential A is pro- 
portional to 1/R, which means that the strength of "ge- 
ometrical" frustration is proportional to the curvature. 

Entry of the first dislocation 

Here we estimate the free energy of a single disloca- 
tion and study the condition for the creation of the first 
dislocation in the system as the strength of frustration is 
increased. We should note that the "penetration depth" 
A is infinite in the present system. In a A = oo (extremely 
type-II) superconductor, as in the case of rotating super- 
fluid , th e lower critical field H c \ depends on the system 
size |21|. We encounter a similar situation here. 

We assume that the dislocation is located at the center 
of the system (x = 0, y = 0) and, instead of treating the 
boundary condition seriously we employ a rather simple 
cutoff procedure [2^|. The phase modulation due to the 
vortex is put as V# = — z x r/|r| 2 , where z is the unit 
vector in z-direction and r = [x,y). Substituting this 
into Eq. |J5J, we obtain the dislocation free energy — 

Knlog^- — 2(/*r W 2 . Here the divergence at y = is 

cutoff at \y\ = d y — d y /"f (d y is the lattice constant in 
j/-direction) . The condition for the dislocation creation 
(namely < 0) is, then, obtained as 
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Many dislocations and destruction of elasticity 

Now we study the strongly frustrated regime where 
many dislocations are created. In this case, the size of 
the dislocation "core" is important. To study this, we 
discrctize again the y-coordinate in Eq. Q as follows, 
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Because of large anisotropy Ai ^> A3, the core is 
strongly elongated along x-direction although the size in 
y-direction is comparable to the lattice spacing. There- 
fore we extract two crystal planes between which the core 
is located and neglect the coupling to other planes. The 
solution which minimizes F in Eq. can be obtained 
analytically as %+i = — Uj — (d x /~ir) arctan(x/A) (the 
core is locate between the j-th and the (j + l)-th plane). 



The core size in x-direction, A, is given by A = d y /j (see 
Ref. [23] for a similar calculation). If we consider that 
the core size in y-direction to be d y , the core size in the 
scaled system (see Eq. ©) is given by d y /j in both x- 
and y-direction and the core is considered to be isotropic. 

In case of a type-II superconductor, the order is de- 
stroyed by a magnetic field when the distance between 
vortices becomes comparable to the vortex core radius, 
and this field strength is called the upper critical field 
Hc2 — </>o/(2tt£ 2 ), where £ is the coherence length. In 
the present case, noting the correspondences, 



H 



7 



d x R 7 



(12) 



where H is the magnetic field, we obtain the critical 
radius R C 2 = 2ndy/(-fd x ) corresponding to H C 2. If 
R < Rc2 , the frustration due to bending is so strong that 
the crystal is full of dislocations, whose cores are almost 
overlapped. In this case, each plane can be almost free 
from the elastic coupling to neighboring planes. In other 
words, the inter-plane elasticity is destroyed by the dislo- 
cations. If R C 2 < R < R c i, the density of dislocations is 
low enough and the cores are separated from each other. 
In this case, the crystal may retain a finite elasticity. 

Phase diagram of ring crystals 

The R-W phase diagram of ring crystals is depicted in 
Fig. 13 As is simply imagined, R > W/2 is an absolute 
requirement for the formation of a ring crystal. The for- 
bidden area is hatched in the figure. It should also be 
noted that our approach is not reliable in the region too 
close to R — W/2, because of the assumption R ^> W. 

We have divided the whole region into three "phases" (I 
~ III). In the phase I, the crystal is free from disloca- 
tions (the "Meissner phase" in superconducting analog). 
The phase boundary R c \ is reentrant as a function of W 
because of the logarithmic contribution in Eq. (|1U|> . al- 
though the lower part of the curve is meaningless because 
W < s/edy. (We denote the smallest R in the phase I 
by -Riim-) In the phase II, the frustration is not too large 
and the system contains dislocations well separated from 
each other (the "mixed state" ) , and the crystal may still 
hold rigidity. In the phase III, the frustration is large 
and the dislocation cores are overlapped, destroying the 
inter-plane rigidity of the crystal (the "normal state"). 
Possible dislocation pattern in each phase is also depicted 
in Fig. 

The distribution of dislocations are observed directly 
by scanning tunneling microscope etc. Phonon spectrum 
measurement is also a powerful tool to investigate the 
crystal rigidity, since there are drastic changes in elastic- 
ity among I ~ III. This may affect the physical properties 
of the system also. Studies on various ring crystals from 
such point of view may be an interesting topic in the fu- 
ture |2J, [25| . Actually, Tanda et al. have observed that 
the critical temperature of charge density wave transition 
in NbSe3 is significantly modified in ring crystals, whose 
relation to crystal defects may be worth investigating. 
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FIG. 3: The phase diagram of a ring crystal and the disloca- 
tion pattern in each phase (see text). In the crystal pictures, 
dislocation cores are shown by parts of arcs. 

Now we apply our theory to actual systems. First, we 
study multi-wall carbon nanotubes, whose elastic con- 
stants are approximately those of graphite. From Ref. 
[2r| . we read Ai = en ~ 100, A 2 = c 33 ~ 4, A 3 = 
c 44 ^ 0.02 (dyn/cm 2 ), and A3/A1 = 7 2 ~ 10~ 4 . Using 
d x = 0.25 nm and d y = 0.34 nm, we obtain R C 2 = 46.2 
nm and Ru m = 395 nm. When R < R&, the crystal 
is full of dislocations, and the elastic coupling between 



graphite sheets is negligible. This is consistent with the 
fact that thin nanotubes consist of graphite sheets with 
different chiralities and the inter-plane coupling becomes 
significant only in thicker tubes (R > 75nm) |27| . How- 
ever an accurate treatment of different chiralities based 
on dislocation theory is a future problem. 

In case of transition metal trichargogenides, we study 
orthorhombic TaS 3 , for which the order of A1/A2 is es- 
timated from Young modulus E and shear modulus G 
in Ref. jH |U as G/E = A3/A1 = 7 2 ~ 10~ 2 . Using 
d x = 0.3 nm and d y = 2 nm, we obtain, R c i — 133.3 
nm and Ru m = 1138 nm. Although the phase I may 
be difficult to realize (e.g., R ~ 30/im and W — 40nm), 
crossover from II to III may be experimentally accessible. 

It is also interesting to investigate tubes with R ~ W 
by numerically minimizing the free energy, Eq. or JfjJ) , 
instead of using superconducting analog. Such study is 
now prepared. 

In summary, we have clarified the distribution of dis- 
locations, in the ring-(or cylinder-)shaped crystals with 
various radiuses and thicknesses. Our results may be 
useful for the studies of the physical properties of "topo- 
logical matters" , and for their future applications to elec- 
trical and mechanical devices. 

MH and HE were financially supported by a Grant- 
In- Aid for Scientific Research from the Ministry of Ed- 
ucation, Science, Sports and Culture, Japan. KK was 
financially supported by the Sumitomo Foundation. 



[1] N. W. Ashcroft and N. D. Mermin, Solid State Physics 

(Saunders college, 1976). [16 
[2] S. Iijima, J. Phys. Chem. p. 3466 (1987). 

[3] R. F. Curl and R. E. Smally, Science 242, 1017 (1988). [17 
[4] H. W. Kroto, Science 242, 1139 (1988). 
[5] S. Iijima, Nature 354, 56 (1991). [18 
[6] R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical 

Properties of Carbon Nanotubes (Imperial College Press, [19 

1998). 

[7] M. Hayashi, Phy. Lett. A 342, 237 (2005). [20 
[8] S. Tanda, H. Kawamoto, M. Shiobara, Y. Okajima, and [21 

K. Yamaya, J. Phys. IV France 9, 379 (1999). 
[9] S. Tanda, H. Kawamoto, M. Shiobara, Y. Sakai, S. Ya- [22 
suzuka, Y. Okajima, and K. Yamaya, Physica B 284- 
288, 1657 (2000). [23 
[10] Y. Okajima, H. Kawamoto, M. Shiobara, K. Matsuda, [24 
S. Tanda, and K. Yamaya, Physica B 284-288, 1659 
(2000). [25 
[11] S. Tanda, T. Tsuneta, Y. Okajima, K. Inagaki, K. Ya- [26 

maya, and N. Hatakenaka, Nature 417, 397 (2002). 
[12] T. Tsuneta and S. Tanda, J. Crystal Growth 264, 223 [27 
(2004). 

[13] G. Seifert, H. Terrones, M. Terrones, G. Jungnickel, and 

T. Frauenheim, Phys. Rev. Lett. 85, 146 (2000). [28 

[14] G. Griiner, Density Waves in Solids (Addison- Wesley, [29 
1994). 

[15] M. Hayashi and H. Ebisawa, J. Phys. Soc. Jpn. 70, 3495 



(2002). 

M. Hayashi, H. Ebisawa, and K. Kuboki, Phys. Rev. B 
72, 024505 (2005). 

K. Yakubo, Y. Avishai, and D. Cohen, Phys. Rev. B 67, 
125319 (2003). 

K. Wakabayashi and K. Harigaya, J. Phys. Soc. Jpn. 72, 
998 (2003). 

L. D. Landau and E. M. Lifshitz, Theory of Elasticity 
(Pergamon, New York, 1986). 

M. Tinkham, Superconductivity (McGraw-Hill, 1996). 
W. F. Vinen, chapt. 20 of Superconductivity, R. D. Parks 
ed., (Marcel Dekker, 1969). 

M. Hayashi and H. Yoshioka, Phys. Rev. Lett. 77, 3403 
(1996). 

M. Hayashi, Physica B 324, 82 (2002). 

K. Shimatake, Y. Toda, and S. Tanda, cond- 

mat/0511328. 

T. Nogawa and K. Nemoto, cond-mat/0511326. 

O. L. Blakslee, D. G. Proctor, E. J. Seldin, G. B. Spence, 

and T. Weng, J. Appl. Phys. 41, 3373 (1970). 

Y. Maniwa, R. Fujiwara, H. Kira, H. Tou, E. Nishibori, 

M. Takata, M. Sakata, A. Fujiwara, X. Zhao, S. Iijima, 

et al., Phys. Rev. B 64, 073105 (2001). 

X. D. Xiang and J. W. Brill, Phy. Rev. B 36, 2969 (1987). 

K. Yamaya and G. Oomi, J. Phys. Soc. Jpn. 51, 3512 

(1982). 



